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September 24, 2009 

ABSTRACT 

It has been hypothesized that large fluxes of neutrinos may be created in 
astrophysical "cosmic accelerators." The primary background for a search for 
astrophysical neutrinos comes from atmospheric neutrinos, which do not exhibit 
the pointlike directional clustering that characterizes a distant astrophysical sig- 
nal. We perform a search for neutrino point sources using the upward-going muon 
data from three phases of operation (SK-I, SK-II, and SK-III) spanning 2623 days 
of live time taken from April 1, 1996 to August 11, 2007. The search looks 
for signals from suspected galactic and extragalactic sources, transient sources, 
and unexpected sources. We find interesting signatures from two objects — RX 
J1713. 7-3946 (97.5% CL) and GRB 991004D (95.3% CL)— but the signatures 
lack compelling statistical significance given trial factors. We set limits on the 
flux and fluence of neutrino point sources above energies of 1.6 GeV. 

Subject headings: neutrino astronomy, Super-Kamiokande, microquasars, pleri- 
ons, gamma-ray bursts, supernova remnants, active galactic nuclei 



1. Introduction 

The most plausible explanation of ultra-high-energy (UHE) cosmic rays is the so-called 
cosmic accelerator model, which posits that protons are accelerated by electromagnetic en- 
ergy associated with objects with extended magnetic fields. It has been pointed out that 
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cosmic accelerators, through pion production, may be capable of producing copious high- 
energy neutrinos (in excess of 1 G e V) at fluxes detectable by current or planned experi- 
ments (jBahcall fc Waxmannll200ll ). 



Many models invoke the physics of shocks, which pr edicts a power law spectrum that 
falls off like d<& v / dE oc E -1 with a spectral index of 7 ~ 2 ([Learned fc Mannheimll2000l ). The 
cross section for neutrino-nucleon scattering, however, depends roughly linearly on energy, 
and the effective detector volume increases with energy since high-energy muons can travel 
further in rock than lower-energy ones. For the most likely range of spectral indices (7 = 
2 ~ 3), the vast majority of astrophysical neutrino events at Super-Kamiokande fall into the 
category of upward-going muons described in Section [2j The spectrum of neutrino-induced 
muon events at Super-Kamiokande is shown in Figure [TJ 



2. Upward-Going Muons at Super-Kamiokande 

Super-Kamiokande is a water Cherenkov detector located in Mt. Ikeno in central Japan 
under 2700 meters water-equivalent rock overburden. The cylindrically shaped detector 
contains 50ktons of ultrapurified water. It is divided by a steel structure into an inner 
detector (ID) equipped with 11146 50 cm PMTs aimed inward and an optically separated 
outer detector (OD) instrumented with 1885 20 cm PMTs aimed outward and equipped with 
wavelength-shifting plastic plates. The primary function of the OD is to identify charged 
particles that enter or exit the ID most of which are cosmic rays. Within the ID we define a 
22.5 kton fiducial volume, within which detector response is expected to be uniform. Further 
details regarding the Super-Kamioka nde detector design, op er ation, calibrations, and data 



reduction can be found in References Fukuda. S. et al. 



(12003 ): lAshie. Y.. et all (120051 ) . 



Neutrino events with total deposited energy above ~ 100 MeV are overwhelmingly due 
to atmospheric neutrinos from the decay of pions created by cosmic rays in the upper at- 
mosphere. To good approximation the Earth is transparent to neutrinos up to energies of 
order of 100 TeV. (The neutrino interactions observed by Super-Kamiokande represent a 
tiny fraction of the 1 GeV — 100 TeV neutrinos passing through the Earth.) 

Super-Kamiokande classifies events according to their energy and/or topology (see Fig- 
ure [2]). Fully-contained (FC) neutrino events are those where interaction products are ob- 
served in the ID, with no significant correlated activity in the OD. Partially-contained (PC) 
events are those where some interaction products exit the ID. Upward-going muon (upmu) 
events occur when a penetrating particle traveling in the upward direction enters and either 
stops in or passes through the detector. Upmu events are attributed to muons produced 
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by neutrino interactions in the rock surrounding the detector. A cut on reconstructed path 
length ensures that upmus have energies greater than E™ n = 1.6 GeV. FC, PC, and upmu 
events are associated with successively higher energy samples of neutrino interactions, rang- 
ing from 200 MeV for the lowest energy FC events to above 1 TeV for the highest energy 
upmus. We may thus think of each event type as a coarse energy bin so long as it is under- 
stood that the spectra for each event type overlap significantly. 

Upmus are subdivided into stopping upmus, which stop inside the detector, and through- 
going, which pass straight through. The through-going upmu sample is further divided into 
showering and non-shower ing upmus with chara cteristic neutrino parent energies of 800 GeV 



and 100 GeV respectively (IDesai. S. et al.ll2007l ). For point sources with spectral indices near 
7 = 2, neutrinos are overwhelmingly likely to be observed in the through-going muon channel, 
and so we henceforth focus our attention on this dataset. The upmu dataset includes a small 
contribution (< 1%) from nearly-horizontal cosmic rays, which constitute an additional 



results making use of the upmu data include ( 


Desai. S. et al.l 


2007; 


Fukuda, S., et al. 


2002; 


Thrane. E.. et al. 


2009; 


Swanson. M. E. C. et al. 


2006: 


Abe. K.. et al. 


2006). We address 



nearly-horizontal cosmic rays again in Section [3J 

We utilize data from three stages of operation summarized in Table [H The SK-II stage 
was characterized by diminished phototube coverage, (about half the value for SK-I and 
SK-III). A comparison of sola r-neutrino data from SK-I and SK-II has found no systematic 



differ ences in the two stages (jCravens. J. P.. et al. (The Super-Kamiokande Collaboration) 



20081 ). The reduced photocoverage during SK-II caused the upmu angular resolution to 
worsen marginally: 1.46° for SK-II in comparison with 1.05° for SK-III. (Angular resolution 
is defined such that 68% of MC events have an angular separation between the true and 
reconstructed muon direction that is less than the resolution.) We take this difference into 
account in Section [6] in the construction of our search algorithm. SK-I and SK-III are similar 
in that they share the same number of PMTs, but they differ in that SK-III has additional 
reflective material called Tyvek® in the OD to optically separate the caps from the barrel 
to aid in event reconstruction. 



3. Likelihood Algorithm 



In this section we develop a likelihood algorithm to maximize th e discovery potentia 



We find it to be more sensitive by a factor of two than a previous search (IDesai. S. et al.112 007') 



wher ein we counted the number of showering events in a 4° search cone, (see also (lAbe. K.. et al 
20061 )). A similar algorithm has recently been presented as a means of searching for point 
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dO^/dE (arbitrary units) vs. log(E) 




2 4 6 



Log(E[GeV]) 



Fig. 1. — The spectra of neutrino-induced upmus for a source with spectral indices 7 = 2 
(dashed) and 7 = 3 (solid) and a minimum upmu energy of E™ n =1.6 GeV. The normalized 
spectra cross in the neighborhood of ~ 200 GeV. 




Fig. 2. — The spectra of different ev ent types at Super-Ka miokande obtained using 100 years 
of atmospheric neutrino MC (from (IDesai. S. et al.l 120071 )). 
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sources with IceCube data (IBraun et al.ll2008l ) 



To begin we construct an 8° search cone around the test direction, which is big enough 
to include over 99% of the point spread function (PSF). We posit that any event, i, falling in 
the search cone has a probability a of being point-source signal, S, and a probability (1 — a) 
of being atmospheric background, B, as in Equation [TJ 

P l = aS(9 l ...) + (l-a)B(9 l ...) (1) 

S(9i...) is the probability density distribution for point-source neutrinos. It depends, among 
other things, on the angular separation, 9, from the search direction. (The ellipses indicate 
that S and B depend on additional variables, which we shall enumerate in section 0) 

Smaller values of 9 are more likely than large ones for the signal distribution. B(9{...) 
is the probability density distribution for the background, which consists of atmospheric 
neutrinos plus a small contribution (< 1%) from nearly- horizontal cosmic ray muons that 
can masquerade as neutrino-induced upmus. Larger values of 9 are more likely than smaller 
ones for the background distribution since there is more phase space near the edge of the 
cone than at the center. 

If there are N events in the search cone, we can define a likelihood function, C, as the 
product over events, i, of each P(a\6i...) as in Equation [2J 

N 

C(a) = HP(a\9i...) (2) 

i 

Now we can interpret a as the fraction of events in the search cone due to signal as in 
Equation [31 when a = there is no signal, when a = 1 the signal is maximal. 

(a) = N s /(N s + N B ) (3) 

Here N$ and Nb are respectively the number of events in the cone due to signal and back- 
ground. 



Table 1. A summary of the dataset for this analysis. (SK-III continued through 2008, but 
the data used here are limited to the dates indicated.) 



phase 


dates 


events 


SK-I 


Apr 1, 1996 - July 19, 2001 


1,879 


SK-II 


Jan 17, 2003 - Oct 5, 2005 


888 


SK-III 


Aug 4, 2006 - Aug 11, 2007 


367 
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As it stands, the likelihood function depends only on the properties of each event (such 
as 8) and not on the excess/deficit of events in the search cone in comparison to the expected 
number of events given a, which we denote Nm, and which is an additional measure of signal 
strength. N(a) is the expected number of events given a whereas N is the number of observed 
events. 

To add this information to the likelihood function, we define a generalized likelihood 
function, denoted £, which is the product of £ with a Poisson probability distribution on a, 
as in Equation |H 

-2(a) = -Ppoisson£(a) = e~ N( - a) — ^ JJ P( a |^.„) (4) 

i=l 

Noting that Nt a \ = N$ + Nb, and recalling Equation [31 we see that Nt a ) = Nb/(1 — a). 
Nb, the expected number of background events, can be determined by MC, and so we can 
eliminate all unknowns from Equation S] in order to obtain Equation 

Having constructed the generalized likelihood function, we are now able to define a 
procedure for assessing signal strength. Given a search direction from which we can define 
{Oi}, we vary a in order to find the best fit value, op, which maximizes the generalized 
likelihood. Finally, we define a likelihood ratio, A, as in Equation [HI which serves as the 
ultimate indicator of signal strength. 

A = 2 log [£{a F )/£(a = 0)] (6) 

As A gets larger, the probability of observing A due to fluctuations in the atmospheric 
background, P{> A), becomes smaller. P{> A) is determined numerically using atmospheric 
MC. 



4. Sensitivity 

In one test described in Section [71 we measure A at regularly spaced 0.5° intervals over 
the entire sky visible to Super-Kamiokande. First we determine the largest signal at any 
point in the sky, A max , and then we ascertain if it is large enough to constitute a statistically 
significant signal using the probability distribution, P(A max ), which is determined with MC. 
Studying P(A max ) (see Figure [3]) we find that the detection of a source at dec = —15° (in 
the middle of the range of declinations visible by Super-Kamiokande) at 90% (99%) CL 
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with an efficiency of 50% requires a signal of AMax = 30.2 (34.7) respectively. We make a 
distinction between the false alarm rate, which determines the threshold for A 90 % (A 99 %), 
and the efficiency, which tells us the probability that a true signal will exceed the threshold. 
These values of A max correspond to 10 (11) upmus out of 3134 total events. Assuming a 
spectral index of 7 = 2, ten upmus corresponds to a neutrino flux of ~ 3 x 10~ 7 cm _2 s _1 
above 1.6 GeV. 

The sensitivity calculated here describes the flux needed to observe a statistically sig- 
nificant signal by looking everywhere in the sky at regularly spaced 0.5° intervals. Such a 
search does not identify any a priori candidate sources, so we call it a "tabula rasa search." 
It is important to note that the sensitivity of a search for a single a priori candidate is sub- 
stantially better than the tabula rasa search because the latter samples the entire visible sky 
where the presence of fluctuations is sure to push the detection threshold higher. 



A previous search for astrophysical neutrinos (IDesai. S. et al.ll2007l ) used the showering 
muon dataset to define search directions. A 4° cone was drawn around each search direction 
and the number of events inside the cone was compared to the number expected from atmo- 
spheric background. The 4° "hard cone search" did not incorporate information about the 
angular separation of the events from the search direction and it did not make use of the 
non-showering dataset. We find that the algorithm presented here requires half the signal of 
the hard cone test to obtain a signal at 90% CL all other things equal. 



5. Upper Limit Calculation 

We also use the generalized likelihood function to calculate upper limits on the number 
of point-source neutrino-induced upmus as in Equation 

ago=ap+8 i>1 

£(«)/ / £(a) = 90% (7) 



' ap— 8 JO 

The quantity ago = (at + 5) represents the upper limit at 90% CL on the fraction of events 
in the search cone due to signal. (In the not infrequent event that (ap — 5) is unphysically 
less than zero, we modify Equation [7] so that the lower limit is 0.) We thereby obtain an 
upper limit on the E v > 1.6 GeV point-source neutrino-induced upmu flux as in Equation [HI 

"90 N , a s 



M An- + 

^J-effective ^ 



effective '-exposure 



Here v4 e g- e ctive is the effective area of the detector and ^exposure is the exposure time. (Both 
^effective and Exposure depend implicitly on the search direction since they each vary with 
zenith angle, z.) We defer discussion of the neutrino flux calculation until Section [91 
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6. Likelihood Variables 

Having introduced the formalism of the likelihood function, we are ready to give a 
full account of the variables on which it depends. In addition to angular separation, 9, the 
likelihood function depends on event type, denoted n, which can take on values of "showering" 
and "non- showering." The atmospheric neutrino flux falls with a steep d<& u /dE oc E~ 3 7 
spectrum in the relevant energy range, whereas point-source neutrinos are thought to fall 
off like d$ v /dE oc E~ 2 . Since showering muons come from a higher energy parent neutrino 
population than non-showering muons, showering muons provide stronger evidence of point- 
source neutrinos than non-showering muons. 

The likelihood function also depends on n since the PSF is different for showering and 
non-showering events. (Showering events, which produce excess light through pair produc- 
tion, bremsstrahlung, and photonuclear effects, are on average slightly harder to reconstruct 
than non-showering events, which deposit energy primarily through ionization.) An example 
of the PSF for showering and non-showering muons is depicted in Figure HI 

The likelihood function also depends on detector geometry, denoted m, which can take 
on values of "SK-I," "SK-II," or "SK-III." The PSF and showering algorithm efficiency both 
vary with m. The likelihood additionally depends on each event's zenith angle, z, as well as 
the equatorial coordinates of the search direction: right ascension, ra, and declination, dec. 
Dependence on z arises from modest differences in the PSF at different zenith angles. (The 
fitter performs best for slightly upward-going going muons and worst for straight upward- 
going muons due to detector geometry.) 

Dependence on (ra,dec) enters through the number of expected background events in 
the cone, Nb- Atmospheric neutrinos have a well-studi ed zenith angle dependence due 



to the "secant theta effect" (ILearned fc Mannheim! 120001 ) from mesons interacting in the 



atmosphere plus additional effects from oscillations and the detector shape. Each equatorial 
coordinate (ra,dec) is associated with a different locus of zenith angles and so the zenith 
angle distribution of atmospheric neutrinos manifests itself as a dependence on equatorial 
coordinates. Thus, each measurement of signal strength is characterized by the following 
variables: 



N 



Y\_{9i,ni\mi, Zi, th, dec} (9) 

i 

In order to determine the signal and background distributions we use MC. The back- 



ground MC begins with atmospheric upmu vectors gener ated with NEUT (iHayatol |2002| ) 



and using the Honda model of atmospheric neutrino flux (IHonda et al.lll995l ). These upmu 



vectors become the input for a GEANT 3 based detector simulation. The same MC has been 
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Fig. 3. — A normalized histogram of A max . The thresholds for 90%/99% CL detection are 
respectively indicated with solid red and dashed green vertical lines. Search directions are 
spaced at 0.5° intervals. 



0.45 r 

0.4 r ... T 




9 (deg) 

Fig. 4. — The normalized point spread function for showering (red squares, dashed line) and 
non-showering (black circles, solid line) muons includes contributions from scattering as well 
as detector resolution. Here we assume a spectral index of 7 = 2. The vertical error bars 
are too small to see. 
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described in previous Super-Kamiokande publications, e.g. lAshie. Y.. et al.l (120051 ). The out- 
put of this simulation is then processed by the upmu reduction as if it were real data. We 
account for neutrino oscillations by applying a weighting factor proportional to the two-flavor 
survival probability, assuming maximal mixing and Am^ = 2.5 x 10~ 3 eV 2 . The procedure 
for accounting for contamination from nearly-horizontal co smic rays using variat ions in the 
rock overburden is detailed elsewhere, (see, e.g., Reference iDesai. S. et al.l (120071 ).) 

To generate signal MC we first assume that the source is characterized by a d§ u /dE oc 
E~"< power law with a spectral index of 7 = 2. We obtain a spectrum for point source 
upmus by multiplying d<§> v jdE by P(E U , E™ m ) — the probability that a neutrino with energy 
E v creates a muon with energy greater than E™ m = 1.6 GeV — and then integrating with 
respect to E v as in Equation [TSJ (We explore this equation in greater detail in Section 
in the context of the neutrino flux calculation.) Each upmu vector is assigned a scattering 
angle according to its energy before it is fed into the detector simulation, which is followed 
by the upmu reduction. 



Four Searches 



We conduct four searches for neutrino point sources. The first search, mentioned in 
Section HJ is the tabula rasa search. For this search we make no assumptions about a priori 
suspected sources. We measure signal strength A at regular 0.5° intervals over the entire sky 
visible to Super-Kamiokande and record the largest observed value, A max , which we compare 
with MC to determine if a statistically significant signal is present. This 0.5° spacing is 
chosen to be smaller than the width of the point spread function (a little over a degree), but 
large enough to allow relatively fast computation. The second search is one for suspected 
candidates. In this search we consider sixteen objects identified in various publications 



as plausible bright sources of astrophysical neutrinos ( 


Zhang et al. 


2003; 


Guetta & Amato 


2002 


Huang & Po 


il 


2008 


Aiello. S.. (The NE 


M 


O Collaboration)! 


20 


A 


Link & Bureio 


2005 


Kistler & Beacom 


2006; 


3ednarek et al. 


2005 




Gabici & Aharonian 


2007 


). Candidates in- 



elude magnetars (young pulsars with strong magnetic fields and periods on the order or 
seconds); plerions (aka pulsar wind nebulae, nebulae with an embedded pulsar); supernova 
remnants (SNR, supernovae Shockwaves colliding with surrounding gas); and microquasars 
(nearby black holes devouring a neighboring star). A summary of the candidates tested is 
provided in Table [2J 

In the third search, we test for a correlation between upmus and 27 ult ra-high- energy cos 



mic rays observed by the Auger experiment, which may be li nked to AGN ([Abraham. J et al. 



20071 ; iHague. J. D. for the Pierre Auger Collaboration!l2009l ). We add the signal strength, A 
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from each direction to create a summed signal, As, the significance of which we evaluate 
with MC. 

There are many ways to search for a correlation of neutrinos with AGN, but this method 
has several advantages. There are more than 300 known AGN within 75 Mpc of Earth. If 
we sum the signal from a great many cataloged AGN, we will effectively sample the entire 
sky, thereby diluting any signal. It is not clear how to pick the most neutrino-bright AGN, 
or how many to include in a summed signal analysis. 

The 27 UHE events observed by Auger, however, constitute a conveniently small dataset. 
Also, it is plausible that the brightest cosmic ray objects are also the brightest in neutrinos. 
Since the Auger experiment is a southern hemisphere experiment looking at downward-going 
cosmic rays, and Super-Kamiokande is a northern hemisphere experiment looking at upward- 
going muons, all 27 UHE events in the Auger dataset point to equatorial coordinates that are 
also visible to Super-Kamiokande. UHE cosmic rays are deflected by electr omagnetic fields , 
but for a reasonable choice of parameters this deflection is small, ~ 1° (ILee et al.l 119951 ). 
Finally, a test for a correlation of upmus with UHE cosmic rays is simple and transparent. 

The fourth search is a se arch for u pmus coincident with appr oximately 2 200 GRBs in the 

BAT SE JPaciesas et al.lll999h . HETE flHETE Mission. The!l2006f). and Swift rtSwift Collaboration. The 



20081 ) catalogs. Following a previously established procedure (IFukuda. S.. et al 



20021) 



we 

employ a ±1, 000 s window centered on the beginning of photonic observations of each GRB. 
This window size allows for any reasonable delay, positive or negative, between neutrino 
emission and photonic emission, while still providing a very small likelihood of random co- 
incidences due to atmospheric neutrinos. Since the upmu rate at Super-Kamiokande is only 
about 1.6 day -1 , the timing cut alone provides a powerful filter to remove atmospheric back- 
ground. If there is a coincidence in time, we next determine if the event falls within the 8° 
search cone. If so, we evaluate the signal strength, A, of the coincident event. The probability 
that the event is due to atmospheric background is given by the product of the probability 
of a coincidence in time, with the probability of a coincidence in space, with the probability 
of measuring a signal strength > A: 



PbG — -ftime -Pspace P{> A) 

As always, we evaluate P(> A) with MC. 



(10) 
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8. Results 

8.1. Tabula Rasa Search 

The tabula rasa search yielded a signal of A max = 19.1, which is well below the detection 
threshold. In fact, this signal is very small given the distribution of A max in Figure [31 which 
prompted us to rerun the algorithm using "boot-strapped" values for local sidereal time. 
"Boot-strapping" in this context is the practice of assigning random values of local sidereal 
time to each event. This has the effect of scrambling the equatorial coordinates associated 
with each event while preserving other properties such as z and n. The second boot-strap 
run yielded A max = 24.9, which, while still below the the detection threshold, is more typical, 
giving us confidence that the small value of A max = 19.1 is a reasonable fluctuation. Since 
no signal was observed, we set upper limits on the flux of point-source neutrinos. 

Upper limits on the flux of point source upmus/neutrinos as a function of declination are 
presented in Figures [5] and [6] respectively. The limits worsen at higher values of declination 
due to reduced exposure time. In Figure [7] we plot the neutrino flux li mits obtained in 



this st udy (black) alongsid e an estimate of those obtained by AMANDA lAbbasi. R. et al. 



(12009( 1 (red) and MACRO lAmbrosio.. et al. (The MACRO Collaboration)! (J200l|) (dotted 



blue). Also included in this plot is the expected detection threshold for this study (dotted 
green) for the case of a 50% detection efficiency. We provide a sky map of the likelihood 
ratio A in Figure [8j in Figure [9] we provide a sky map of p(A) measured in sigma, and in 
Figure [TU] we present the through-going muon dataset in the form of a sky map. 



8.2. Suspected Candidates 



Of the 16 candidates listed in Table [21 we found a statistically interesting signature 
from one, the SNR RX J1713. 7-3946, which yielded a signal of A = 8.5. Of the 26 events 
in the search cone centered on SNR RX J1713. 7-3946, approximately 5 are attributable 
to signal. RX J1713 .7-3946 is a one of the few known TeV shell-type SNR. It is approx- 
imat ely 1 kpc away (IKoyama et al.l 119971 ; iFukui et al.l 120031 ) and approximately 1600 years 
old (I Wang k, Chenlll997l ). Using MC, we find the probability of an accidental signal of this 
magnitude to be 0.16%, (see Figure [TTJ) Taking into account the fact that we performed 
sixteen tests, the probability of a chance occurrence becomes 2.5%. Limits on the neutrino 
flux from each candidate (assuming a 7 = 2 power law) are recorded in Table [2J 
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Fig. 5. — 90% CL limits on the E v > 1.6 GeV flux of upward-going muons from point- 
source neutrinos as a function of declination. Error bars (too small to see) reflect statistical 
uncertainty created by averaging over ra. 
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Fig. 6. — 90% CL limits on the E v > 1.6 GeV flux of point-source neutrinos as a function of 
declination. Error bars (too small to see) reflect statistical uncertainty created by averaging 
over ra. 
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Fig. 7. — A comparison of the neutrino flux limits (at 90% CL for E„ > 1.6 GeV) 



obtai ned here (black data points) with those obtained by AMANDA (lAbbasi. R. et al. 
20091 ) for E„ > 1.9 TeV (approximated by the dotted red line), and limits by 



MACRO jAmbrosio.. et al. (The MACRO Collaboration)! hooil ) for E u >~ 1 GeV (approx- 
imated by the dashed blue). The dotted green line represents the approximate flux required 
to produce a signal at > 90% CL with a 50% detection efficiency. 




Fig. 8. — A sky map of A in equatorial coordinates. The maximum value (A max = 19.1) is 
below the detection threshold. 
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Fig. 10. — A sky map in equatorial coordinates depicting each through-going muon event as 
a dot. Red asterisks are showering muons, blue dots are non- showering. 
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Fig. 11. — A normalized histogram of A generated for atmospheric (background) events 
falling inside a search cone centered on RX J1713. 7-3946. The coincident event was measured 
to have a likelihood ratio of A = 8.5 (marked with a red line), which has a 0.2% chance of 
being background. 

Table 2. A list of 16 candidates considered to be plausible neutrino point sources with 
associated limits (at 90% CL) on E u > 1.6 GeV neutrino flux. The abbreviation "mag" is 
for "magnetar," "pier" is for "plerion," and "MQ" is for "micro quasar." Estimating the 
distance to these sources can be difficult and inexact as distance estimates often depend on 
measurements of what are assumed to be neighboring objects. The objects in this table are 
associated with estimated distances that range from 0.3 — 50kpc. 



source 


type 


(ra,dec) 


<J>» % (cm-Ss- 1 


) 


SGR 1900+14 


mag 


(286.8° 


+9.3°) 


1.12 + 0.12 x 10 


-7 


SGR 0526-66 


mag 


(81.5°,- 


-66.0°) 


1.15 + 0.13 x 10 


-7 


IE 1048.1-5937 


mag 


(162.5° 


-59.9°) 


6.71 + 0.74 x 10 


-8 


SGR 1806-20 


mag 


(272.2° 


-20.4°) 


1.67 + 0.18 x 10 


-7 


Crab 


pier 


(83.6°, +22.0°) 


1.66 + 0.18 x 10 


-7 


Vela X 


pier 


(128.5° 


-45.8°) 


6.87 + 0.76 x 10 


-8 


G343.1-2.3 


pier 


(257.0° 


-44.3°) 


6.81 + 0.75 x 10 


-8 


MSH15-52 


pier 


(228.5° 


-59.1°) 


1.12 + 0.12 x 10 


-7 


RX J1713. 7-3946 


SNR 


(258.4° 


-39.8°) 


2.67 + 0.29 x 10 


-7 


Vela Jr. 


SNR 


(133.2° 


-46.3°) 


9.16 ± 1.0 x 10" 


8 


MGRO J2019+37 


SNR 


(305.2° 


+36.8°) 


2.46 + 0.27 x 10 


-7 


SS433 


MQ 


(288.0° 


+5.0°) 


1.16 + 0.13 x 10 


-7 


GX339-4 


MQ 


(255.7° 


-48.8°) 


5.50 + 0.61 x 10 


-8 


Cygnus X-3 


MQ 


(308.1° 


+40.8°) 


1.32 + 0.15 x 10 


-7 


GRO J1655-40 


MQ 


(253.5° 


-39.8°) 


1.26 + 0.14 x 10 


-7 


XTE J1118+480 


MQ 


(169.5° 


+48.0°) 


1.29 + 0.14 x 10 


-7 
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8.3. Correlation with UHE Cosmic Rays 

We found the summed signal from the directions of 27 UHE cosmic rays possibly linked 
to AGN by the Auger experiment to be As = 10.3, which does not constitute a statistically 
significant signal. We therefore set a limit on the E v > 1.6 GeV flux of neutrinos from Auger 
UHE directions by averaging the limits obtained for each direction, and found: $^ 0% = 
1.06 ±0.12 x lO^cm-V 1 . 



8.4. Correlation with Gamma-Ray Bursts 

Of the 2200 cataloged GRBs considered, 971 met the criteria of occurring during 
Super-Kamiokande live time with a zenith angle no greater than 8° above the horizon. 
The data for this search does not include the recent so-called "naked-eye" GRB 0803 19B, 
which was the brightest GRB observed to date. Due to i ts unusual intensity, a sepa- 



rate search was performed fo r neutrinos from GRB 080318B (IThrane. E.. et al. 1 120091 ) . see 



also flAbbasi. R.. et all 12009 ). 



Of the 971 GRBs visible to Super-Kamiokande, we observed one coincident GRB (desig 
nated 991004D) with the upmu dataset. The upmu occurred 411 s after the G RB and was sep- 



arate d by 3.4°. GRB 991004D was a long GRB with a duration of T90 = 34 s (IFrontera et al. 



20091 ). The total and peak gamma-ray fluence were not atypical nor was the spectral hard- 



ness. Key properties of the coincident GRB are listed in Table El 

There were 3134 upmus considered in this dataset with a live time of 2.3 x 10 8 s, which 
implies the mean time between upmu events is r = 7.3 x 10 4 s. Thus, the Poisson probability 
for detecting one or more background events in the t = 2000 s window is given by: 

P timc = 1 - P( n = 0) = 1 - e~ tjT = 0.0270 (11) 

The search used an 8° half-angle search cone. Taking into account the fact that the 
search cone sometimes overlaps with the insensitive region (z < 0), the effective search cone 
angle is 6.7°. The approximate probability of a coincident background event falling inside 

Table 3. Details of GRB 991004D. 



ID ra dec UT JST=UT+9 



GRB 991004D 210.75° 



-19.53° 13:13:21 22:13:21 
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the cone is given by: 

Ppacc = ficonc/^SK = 0.043/11.37 = 0.00375 (12) 

(Here Q$k is the solid angle of the sky visible to Super-Kamiokande.) Thus, the probability 
of observing one or more background events within the 8° cone is P C oincidence = -Ppace -ftime = 
0.0038-0.027 = 1.03 x 10~ 4 . 

The signal strength for the coincident upmu was determined to be A = 4.6. From a MC 
sample of random background events (see Figure [T2l we determine that the probability of 
observing A > 4.6 with background is Pa = 0.486. 

It follows that the probability for observing a random coincidence with A > 4.6 in a 
single trial is Ptrfai = Pa -^coincidence = 5.00 x 10~ 5 . The probability of observing no coincidence 
with A > 8.5 in one trial is, of course, 1 — P t riai = 0.99995. If we test 971 GRBs, the 
probability of observing no coincidences with A > 4.6 in any of the 987 trials is given by 
(0.99995) 971 = 0.953. Thus, the probability of observing one or more coincidences with 
A > 4.6 during 971 trials is 1 — 0.953 = 0.047. Thus, testing 971 GRBs against the upmu 
dataset, the single observed GRB-upmu coincidence has a 4.7% probability of being due to 
random background. If we had used a smaller ±500 s window instead, the confidence level 
would improve from 95.3% to 97.6%. 



A previous study by Super-Kamiokande (IFukuda. S.. et al.l |2002j) found no evidence 



for a correlation of neutrinos with GRBs and set a E u > 1.6 GeV fluenc e uppe r limit of 



_p90% > o.038cm _2 . Using a 15° search cone, Reference IFukuda. S.. et al.l (120021 ) reported 
a single coincidence — also GRB 991004D — but was unable to infer a statistically significant 
signal. A likely reason for the apparent discrepancy in significance is that the algorithm 
presented here uses a smaller cone in concert with a likelihood function to filter out events 
that contribute to the expected background of the 15° cone search. 

The best limits on GRB fluence currently come from the AMANDA experiment, which 
reports an upper limit of 1.4 x 10~ 5 cm" 2 for neutrino energies between 250 — 10 7 GeV and 



assuming a spectral index of 7 = 2 (lAchterberg. A. et al.ll2008l ). In their search, AMANDA 



employs a timing window defined as the T90 start and end times for the burst during which 
90% of the total background-subtracted counts are observed. (The interval begins when 
5% of the counts have been observed.) Typically T go is on the order of seconds. While 
such a timing window is motivated by the plausible assumptions that GRB neutrinos are 
emitted at the same times as photons, (and that their journey to Earth is unaltered by any 
new physics), it is nonetheless true that such a narrow window would not have allowed for 
detection of the 411 s delayed upmu observed in coincidence with GRB 991004D. The limits 
calculated by AMANDA, therefore, do not apply to models of GRBs (or new physics) that 
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predict delayed or early neutrino arrival times. 



Given AMANDA's limits, along with the fact that there is little in the literature to 
suggest that GRB 991004D was in any way unusual, we interpret the observed signal as a 
background fluctuation and thereby obtain a limit on the average E v > 1.6 GeV fluence of 
neutrinos from GRBs to be: F® 0% = 0.060 ±0.007 cm . This limit is slightly wo rse than the 
one obtained in the previous Super-Kamiokande study (IFukuda. S.. et al.ll2002l ) due to the 
fact that a nearly significant signal necessarily raises the upper limit. Since the new search 
algorithm employed here is more sensitive, the apparent signal from GRB 991004D is bigger, 
and so our limit actually worsens. 



9. Neutrino Flux Limits 



I n order to compute neutrino flux limits from muon flu x limits, we follow (IThrane. E.. et al 



20091 ; ISwanson. M. E. C. et all2006l ; iDesai. S. et al.ll2007l ) and assume a model for the source 
spectrum. It is typically assumed that d$ u /dE oc E~"< with 7 « 2. Here we assume 7 = 2 
and investigate the implications of this assumption in Section [TTJ Given this assumption, the 
flux of neutrino-induced muons is peaked at high energies, between 1 GeV — 1 TeV depending 
on the spectral index (see Figured!) 



Upmu flux is related to the neutrino flux as in Equation 



dE u P(E U} Ef a )3(z, E v 



dE„ 



(13) 



Here $„,(> E™ m ) is the flux of upmus with energies above the minimum upmu energy of 



1.6 GeV. P{E U , E™ m ) is the probability that a neutrino with energy E v creates 
a muon with energy greater than E™ in and S(z, E u ) is the Earth's shadow factor. We use 
cross sections from the GRV94 parto n distribution fun c tion (Gliick et al.lll995l) and the muon 



range is determined using Reference iLipari fc Stanevl (1199 IT ) 



Reno 



P(E, 

Here da C c I 'dE^E,^ E, 



771m] 
"1 fi 



dE, 



da C c 



(E^, E v ) R{E I 



771mm \ 

A" fi ) 



(14) 



iq dE/j 

, is the charged current component of the neutrino-nucleon cross 

section; (neutral current interactions do not produce muons). R{E^, E™ m ), meanwhile, is 
the average range in rock for a muon with an initial energy of E^ and a final energy greater 
than E™ in . 



S(Z, E v ) = e ~^)a(E,)N A 



(15) 
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Here l co i{z) is the Earth's zenith angle-dependent column depth, and Na is Avogadro's 
number scaled by the density of water. The co l umn depth is calculated using the "Prelimi- 
nary Earth Model" in Reference iGandhi et al.l (119961 ). 



10. Systematic Error 

Two sources of systematic error are dominant in our calculation of the neutrino flux. 
One source is due to uncertainty in the neutrino-nucleon cross secti on (used in Equation [T3|) . 



which has an associated uncertainty of ~ 10% at upmu energies (IGandhi et al.lll996l ). To 
ascertain how this affects the neutrino flux uncertainty we calculate the flux two ways, 
$ + = $(er — > 110% a), $_ = $(cr — > 90% a), and estimate the uncertainty as half the 
difference. 

Another systematic uncertainty arises from the fact that we only include charged cur- 
rent interactions in our calculation of the Earth shadow. While neutrinos interacting via 
charged current interactions produce muons and disappear from the incoming beam, neu- 
trinos interacting via neutral current interactions can lose energy but remain in the beam. 
By including only charged current contributions in the Earth shadow, our calculation is too 
optimistic because neutral current interactions reduce the measured flux at the detector. If 
we include both charged and neutral current contributions, however, our calculation is too 
pessimistic since some neutrinos will reach the detector even after undergoing neutral current 
interactions. 

Therefore, we calculate $ two ways: $cc* = — ^cc) an d &nc = — vcc + cnc)- 
The neutrino fluence is estimated as the average of $>cc an d ^cc+nc an d the uncertainty is 
estimated as half the difference. We add these two uncertainties in quadrature and find the 
total uncertainty in <£>„ to be 11%. 



11. Assumptions in the Flux Calculation 

Before concluding, we assess two assumptions made in the calculation of neutrino flux. 
First, we assume a spectral index of 7 = 2, and so it is illuminating to address how the 
flux limits change if we chose a less favorable spectral index. Repeating the neutrino flux 
calculations with a more pessimistic spectral index of 7 = 3, we find that the flux limits are 
raised by approximately three orders of magnitude due to the reduced neutrino-nucleon cross 
section at lower energies. We note that the two spectral indices yield comparable limits for 
the flux of ~ 200 GeV neutrino-induced muons where the normalized spectra intersect (see 
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Figure[T]). Thus, numerical estimates of neutrino flux are extremely sensitive to assumptions 
about the source spectrum. While it is necessary to make an assumption about the nature of 
the source — and the 7 = 2 power law spectrum is certainly a useful standard assumption — it 
should be noted that until we know more about point-source spectra, numerical limits on 
neutrino flux are best suited for comparison between experiments. 

A second assumption implicit in our flux calculations is that point sources emit neutrinos 
and not antineutrinos. The cross section for antineutrino-nucleon interactions is less than 
the cross section for neutrino-nucleon interactions due to the V — A nature of the weak force, 
and this causes the limits for antineutrino flux to be 30% higher than limits obtained for 
neutrinos. 



12. Conclusions 

We have constructed an algorithm for the purpose of detecting neutrino point sources. 
We find that this algorithm is more sensitive by a factor of two in comparison to a previous 
point-source search. Using the algorithm, we performed five tests. We set limits on E v > 
1.6 GeV neutrino flux at 90% CL as low as 6 x 10~ 8 cm _2 s _1 for declinations less than +54°. 
While not competitive with AM ANDA's limits, these limits are among the best for dec < 0°. 
In a test of suspected sources, we looked for signals from 16 predetermined objects classified 
as potentially bright neutrino point sources. Of the 16 objects tested, we found one to have 
a high enough signal to warrant interest. The SNR, RX J1713. 7-3946, was found to possess 
a signature at 97.5% CL. Northern hemisphere neutrino observatories such as ANTARES, 
NEMO, and NESTOR should be able to assess the possibility that the signal is real since 
the object is located at dec = —39.8°. 

Motivated by interest in correlations of AGN with UHE cosmic rays, generated by 
recent Auger experiment results, we looked for a correlation of upmus with UHE cosmic 
rays. We found no evidence for a correlation and set a limit on the average E v > 1.6 GeV 
neutrino flux from Auger UHE event directions: $^ 0% = 1.06 ± 0.12 x 10~ 7 cm -2 s -1 . We 
performed a systematic search for correlations of upmus with GRBs in the BATSE, HETE, 
and Swift catalogs. We found one coincidence (with GRB 991004D), which constitutes 
a signal at 95.3% CL. We set a limit on average GRB E v > 1.6 GeV neutrino fluence, 
= 0.060 ± 0.007 cm- 2 . 
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Fig. 12. — A histogram of A generated for atmospheric (background) events falling inside a 
search cone centered on a GRB coincidence at (ra,dec) = (210.75°, —19.04°). The coincident 
event was measured to have a likelihood ratio of A = 4.6 (marked with a red line), which 
has a 48.6% chance of being background. 



